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1. Introduction 

Taking into account the relevant experimen- 
tal constraints, the CMSSM and NUHMl pre- 
dict that the lightest Higgs boson should have 
couplings similar to those of the Standard Model 
(SM) Higgs boson [iHS] , and that it should weigh 
no more than ^ 130 GeV [IHS]. We recently re- 
ported the results of global frequcntist fits within 
the CMSSM and NUHMl to the first - 1/fb of 
LHC data, also including precision electroweak 
and fiavour measurements and the XENONIOO 
upper limit on elastic spin-independent dark mat- 
ter scattering [3], updating the results of pre- 
vious global fits by ourselves [THH] and oth- 
ers [T51[TS] (see also [17]). The results re- 
ported in [3] included likelihood contours in the 
(mo,?Tii/2), (18-11/3,777.1/2) Sind (Af^, tan/3) planes 
of the CMSSM and NUHMl, as weU as Ax^ func- 
tions for nig, BK{Bs /i^/i^), Mh,MA and 
sparticle production thresholds in e'^e" annihi- 
lation. 

Notable predictions of these global fits included 
Mh = 119.1^2:9 GeV in the CMSSM and Mh = 
nS.Stli GeV in the NUHMl (which should be 
combined with an estimated theory error AM/i ~ 
±1.5 GeV). These two fits are based solely on the 
Higgs- independent searches including the — 
constraint, i.e., they do not rely on the exist- 
ing limits from LEP [18l[19], the Tevatron [20] . 
or the LHC [2T1[22]. These predictions increase 
to Mh = 124.8t?o*5 GeV in the CMSSM and 
126.6t?:^ GeV in the NUHMl if the (5 - 2),, con- 
straint is dropped. 

Subsequently, the ATLAS and CMS Collabo- 
rations have released their official combination of 
the searches for a SM Higgs boson with the first ^ 
1/fb of LHC luminosity at E^m = 7 TeV [23]. Im- 
pressively, the combination excludes a SM Higgs 
boson with a mass between 141 and 476 GeV. 
Most recently, the ATLAS and CMS Collabora- 
tions have presented preliminary updates of their 
results with ^ 5/fb of data [24]. These results 
may be compatible with a SM-like Higgs boson 
around Mh ~ 125 GeV, though CMS also report 
an excess at Mh ~ 119 GeV in the ZZ* channel. 
We recall that, for low values of Mh, the SM elec- 
troweak vacuum would be unstable |25] , decaying 



into a state with Higgs vev > 10^(10^°) GeV if 
Mh = 119(125) GeV, and that a very plausible 
mechanism for stabilizing the vacuum is super- 
symmetry (SUSY) [ig. 

In this paper, we first report the likelihood 
function for an LHC measurement of Mh with 
a nominal value G (115, 130) GeV, incorporat- 
ing the theoretical error ±1.5 GeV and an esti- 
mate ±1 GeV of the possible experimental error. 
In both the CMSSM and NUHMl, this likeli- 
hood function is minimized for Mh ~ 119 GeV 
if (g — 2)^ is included, and is contained within 
the theoretical uncertainty range shown previ- 
ously as a 'red band' [J. We then discuss the con- 
sequences of combining a measurement of Mh — 
125 GeV (assuming that the current excess will 
be confirmed with more integrated luminosity) 
with our previous analysis [3] of constraints on 
the CMSSM and NUHMl including {g - 2)^. 

We find that the best-fit values of mo and 777^/2 
in the CMSSM and NUHMl are moved to sub- 
stantially higher values, especially in the case of 
"71/2- We also update our results on the best- 
fit regions in the (7771/2, tan /3) and (Afyi,tan/3) 
planes, where we find again the substantial in- 
crease in 777i/2, as compared with our pre-LHC 
Mh results. We present the corresponding one- 
dimensional likelihood functions for the gluino 
mass 777g, an average right-handed squark mass 
777q^ , the lighter scalar tau mass, 777^^, as well as 
in the (777^0, cTp^) plane, where 777^0 is the mass 
of the lightest neutralino and is the spin- 
independent dark matter scattering cross section. 
As could be expected, we find larger values of 
777g , 777g^ , 777^0 and 777f ^ than in our pre-LHC Mh 
fit, and smaller values of a^^, though BR(Ss ^• 
fi^lj,^) is little affected. 

Since Mh — 125 GeV is the value that was 
favoured in the CMSSM/NUHMl fits omitting 
the {g — 2)^ constraint [3], we also show some 
results for fits where {g — 2)^ is dropped. In 
this case, we find that preferred regions of the 
(7770,7771/2) planes arc localized at relatively high 
values, corresponding to relatively large sparticle 
masses. Correspondingly, the spin-independent 
dark matter scattering cross section a^^ would 
be relatively small in this case, though again 
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there would be relatively little effect on BR(£?s 

Finally, we show selected results for a hypo- 
thetical measurement of Mh ~ 119 GeV. 

2. Prediction for 

We recall that the independent parameters 
of the CMSSM [17] may be taken as the 
common values of the scalar and fermionic 
supcrsymmetry-breaking masses 'm{),mi/2 at the 
GUT scale, the supposedly universal trilinear soft 
supcrsymmetry-breaking parameter, Aq, and the 
ratio of Higgs v.e.v.'s, tan/3. A study of the dis- 
tribution of Higgs masses in the CMSSM was per- 
formed in [23. Motivated by {g — 2)^ and, to a 
lesser extent, BR(6 — >■ 57), we assume that the 
Higgs mixing parameter /i > 0. In the case of the 
NUHMl [5^, we relax the universality assump- 
tion for the soft supcrsymmetry-breaking contri- 
butions to the two Higgs masses, ra^u^ ~ "'^ffd ' 
allowing m^u^ — rnj^ 7^ uiq. 

In our previous papers [BllTMlBj we constructed 
a global likelihood function that receives con- 
tributions from clcctroweak precision observ- 
ables, B-decay measurements, the XENONIOO 
direct search for dark matter scattering [30] 
and the LHC searches for supcrsymmetric sig- 
nals, calculated within the MasterCode frame- 
work [131 . This incorporates code based on [21] 
as well as SoftSUSY [32], FeynHiggs 2.8.6 [5], 
SuFla [35], Superlso 03], MicrOMEGAs [3S] and 
SSARD [36], using the SUSY Les Houches Ac- 
cord [37]. As before, we use a Markov Chain 
Monte Carlo (MCMC) approach to sample the 
parameter spaces of supcrsymmetric models, and 
the results of this paper are based on the sam- 
ple of TOM CMSSM points and another 125M 
NUHMl points, both extending up to mg, = 
4000 GeV. 

We used in [5] the public results of searches for 
supersymmetric signals using ~ 1/fb of LHC data 
analyzed by the ATLAS and CMS Collaborations 
and ~ 0.3/fb of data analyzed by the LHCb Col- 
laboration. These include searches for jets -I- 
events without leptons by ATLAS [35] and 
CMS [3S], searches for the heavier MSSM Higgs 
bosons, H/A [3TJ[22], and new upper limits on 



BR{Bs M^Ai") from the CMS [3D], LHCb [31] 
and CDF CoUaborations [31]. Our global fre- 
quentist fit [3] yielded regions of the CMSSM and 
NUHMl parameter spaces that are preferred at 
the 68 and 95% CL. 

This was the basis in [3] for the predictions 
Mh = 119.112:9 GeV in the CMSSM and Mh = 
118. 8t^;^ GeV in the NUHMl, if the {g - 2)^ 
constraint is included as calculated using the 
FeynHiggs code which is quoted as having a the- 
oretical error ±1.5 GeV [S]. It is important to 
note that these best-fit values are well above the 
LEP lower limit and below the Tevatron/LHC up- 
per limit on Mh, which played no role in their 
determination. Fig. 12 of [5] displayed the Ax^ 
likelihood functions for the FeynHiggs value of 
Alh in these models as blue lines, with the theo- 
retical error ±1.5 GeV represented by red bands 
in these plots. As already noted, these predic- 
tions increase to Mh = 124.8110*5 GeV in the 
CMSSM and 126.6t";9 GcV in the NUHMl if the 
{g — 2)^ constraint is dropped. The uncertainty 
on the Mh prediction is somewhat asymmetric, 
which is due to the different constraints coming 
into play. At low Mh values, the most impor- 
tant constraint is that due to the LHC, though 
other low-energy constraints also play roles. On 
the other hand, at high values of Mh, it rises log- 
arithmically with the scalar top masses, so in- 
creasing Mh increases exponentially the required 
supcrsymmetry-breaking mass scales, and wors- 
ens the agreement with other low-energy data and 
the CDM constraint. 

Results without a Higgs-boson mass mea- 
surement 

Within the supersymmetric frameworks discussed 
here, a confirmation of the excess reported by AT- 
LAS and CMS [25 and consequently the discov- 
ery of a SM-Hke Higgs boson is expected to be 
possible in the coming year, with a mass in the 
range between 114 and 130 GeV [24] . We assume 
that this measurement will yield a nominal value 
of Mh within this range, with an experimental 
error that we estimate as ±1 GeV. We now esti- 
mate the one-dimensional likelihood function for 
the nominal central value of Mh, which may be 
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written as = M^" + AM^*^ + AM^^'p, where 
M™ denotes the output of FeynHiggs (which 
was plotted in Fig. 12 of [3] for the fits includ- 
ing {g — 2)^), AM^'^ denotes its difference from 
the true value of Mh (the theoretical error es- 
timated as ±1.5 GeV), and AM^^^ denotes the 
experimental error in measuring (estimated 
as ±1 GeV). Here we treat the experimental and 
the theoretical errors as Gaussians, and include 
them as supplementary uncertainties in the fit 
for the nominal central value of Af/j. As a conse- 
quence of including these uncertainties, the Ax^ 
function for the nominal central value of pre- 
sented here differs slightly from the Ax^ function 
for the FeynHiggs estimate M™ shown in Fig. 12 
of 12]. 

We see in Fig. [TlH that the values of Ax^ for the 
nominal value of Mh calculated in the CMSSM 
and NUHMl with the {g — 2)^ constraint and 
including both the theoretical and experimental 
errors lie below the blue lines taken from Fig. 12 
of [3]. This is to be expected, since the calcu- 
lation of the dashed line incorporates additional 
uncertainties. As is also to be expected, in each 
case the calculated Ax^ lies within the previous 
red band. The most likely nominal value of the 
LHC measurement of Mh remains Mh — 119 GeV 
in both the CMSSM and NUHMl. A value of 
Mh ~ 125 GeV is disfavoured in our analysis by 
Ax^ = 2.2 in the CMSSM and by 1.6 in the 
NUHMl if {g — 2)^ is included. For comparison, 
a nominal value of Alh = 114 GeV, correspond- 
ing roughly to the lower limit set by LEP for an 
SM-like Higgs boson HHlIS], has Ax^ = 0.8(1.5). 
On the other hand, if we drop (g — 2)^ there is 
essentially no price to be paid by including a 
measurement of Mh — 125 GeV. 



^The 'theoretically inaccessible' area of > 130 GeV 
could in principle be extended to higher values if one ex- 
tended the scanned ranges of mi/2 "^Oi which are 
both restricted here to below ~ 4 TeV, as discussed above. 
However, due to the logarithmic dependence of M^, one 
would gain only about one GeV even if values up to 10 TeV 
were included. 



3. Implementation of the LHC Constraint 
on Mh 

We now study the possibility that the LHC 
experiments confirm the excess reported around 
125 GeV and indeed discover a SM-like Higgs bo- 
son. Assuming 

Mh = 125 ± l(exp.) ± 1.5(theo.) GeV , (1) 

we incorporate this new constraint using the 'af- 
terburner' approach discussed previously [3] for 
other observables. This value would be favoured 
if {g — 2)^ were dropped from our global CMSSM 
or NUHMl fit [5]. Alternatively, a measurement 
of such a high Mh value could point to the re- 
alization of some different (possibly GUT-based) 
version of the MSSM (see, for instance, [ISj). We 
also mention briefly some implications if Mh — 
119 GeV. 

Comments on the LHC data 

As a preamble to these studies, we first comment 
on the resuhs of the current ATLAS/CMS Higgs 
combination. We recall that the local p-value 
for the background-only hypothesis for the excess 
found in the ATLAS data at Mh ~ 125 GeV is 
p = 1.9 X 10""', while that in the CMS data at 
Mh ~ 125 GeV has p = 5 x lO'^. In addition, 
CMS reports an excess in the ZZ* channel at 
Alh = 119 GeV with similar significance, but this 
is not confirmed by ATLAS. 

In order to assess the global p- value of a poten- 
tial signal, one should take the 'look-elsewhere ef- 
fect' (LEE) into account. This is conventionally 
estimated by adding to the local p- value the quan- 
tity A^x exp(— Z^jj^/2), where N is the number 
of times the observed upper limit on the signal 
crosses over the /i = (t/ctsm = level in the up- 
ward direction, and .^max is the maximal signal 
significance [24]. Accounting for the LEE, AT- 
LAS assess the global value of their excess at 
125 GeV In the range (110, 146) GeV to be 0.6%, 
and CMS assess the significance of their excess at 
125 GeV to be 1.9% in the range (110, 145) GeV. 

On the other hand, as the CMSSM and 
NUHMl naturally require Mh < 130 GeV, the 
LEE factor is strongly reduced in these frame- 



5 




Figure 1. The one- dimensional Ax^ functions for Mh in the CMSSM (left) and the NUHMl (right). 
The solid lines are for fits including all the available data including [g — 2)^ hut excluding the direct 
LEP \18\ \19 \l . Tevatron \20}/ and earlier LHC \21\ \2S ^ constraints on Mh, with a red band indicating the 
estimated theoretical uncertainty in the calculation of Mh of 1.5 GeV. The dashed line shows the Ax^ 
likelihood function for the nominal central value of a hypothetical LHC measurement of Mh, as estimated 
on the basis of the frequentist analysis of J3jl, and allowing for an experimental error of ±1 GcV in the 
measurement of Mh and a theoretical error of ±1.5 GcV in the FeynHiggs calculation of Mh at any given 
point in the parameter space. 



works. 

Since the excess around 125 GeV is common 
to both experiments and has the correct sig- 
nal strength: fi 1 can be interpreted as a 
Higgs signal in either the SM or a supersymmetric 
framework. We focus here on this interpretation, 
commenting subsequently on some imphcations if 
Mh =i 119 GeV. 

What if Mh = 125 GeV? 

We first examine the effects on the global likeli- 
hood functions in various CMSSM and NUHMl 
parameter planes, and then study implications 
for various obscrvablcs of a potential LHC mea- 
surement Mh ~ 125 GeV, see Eq. (P). The 
(TOo,mi/2) planes shown in Fig. [2] are for the 
CMSSM (left) and NUHMl (right) 0. The re- 

^We have omitted from the NUHMl sample displayed here 
and in subsequent figures a grouping of points at large 
S'lid small mo for which different codes yield dis- 
crepant values of the relic density. MicrOMEGAs 1351 and 



gions preferred at the 68% CL are outlined in red, 
and those favoured at the 95% CL are outlined 
in blue. The solid (dotted) lines include (omit) 
the assumed LHC Higgs constraint. The open 
green star denotes the pre- Higgs best- fit point [3] , 
whereas the solid green star indicates the new 
best-fit point incorporating a Higgs-boson mass 
measurement at 125 GeV. 

Since in the CMSSM and NUHMl the radiative 
corrections contributing to the value of Mh are 
sensitive primarily to mi/2 and tan /?, and only to 
a lesser extent to mo (the stop masses, which are 
the most relevant for Mh, depend mostly on mi/2 
due to the RGE running, and only mildly on toq). 



DarkSusy 1441 yield densities within the WMAP range for 
these points, whereas Superlso Relic 1451 and SSARD I36| 
both yield substantially lower densities. The other figures 
shown in this paper are not affected significantly by the 
omission of these points (which have Ax^ > 5), pending 
resolution of this discrepancy. Tests in other regions of 
the NUHMl sample have not revealed significant discrep- 
ancies between these codes. 
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Figure 2. The (TOo,"^l/2) planes in the CMSSM (left) and the NUHMl (right). The Ax^ = 2.30 and 
5.99 contours, commonly interpreted as the boundaries of the 68 and 95% CL regions, are indicated in 
red and blue, respectively, the solid lines including the hypothetical LHC measurement Mh = 125 ± 1 GeV 
and allowing for a theoretical error ±1.5 GeV, and the dotted lines showing the contours found previously 
in ]^ without this Mh constraint. Here the open green stars denote the pre-Higgs best- fit points J^, 
whereas the solid green stars indicate the new best-fit points. 



we expect that the primary effect of imposing the 
Mh constraint should be to affect the preferred 
ranges of mi/2 and tan/3, with a lesser effect on 
the preferred range of mo. This effect is indeed 
seen in both panels of Fig. O We see that the 
68% CL ranges of extend to somewhat larger 
values and with a wider spread than the pre-Higgs 
results, particularly in the NUHMl. However, 
the NUHMl best-fit value of mi/2 remains at a 
relatively low value of ~ 800 GeV, whereas the 
best- fit value of toi/2 in the CMSSM moves to ^ 
1900 GeV. This jump reflects the fiatness of the 
likelihood function for toi/2 between ^ 700 GeV 
and ~ 2 TeV, which is also reflected later in the 
one-dimensional Ax^ functions for some sparticle 
masses 0. 

When we add the hypothetical Mh constraint 
the total the best-fit points increases sub- 

stantially, as seen in Table [1] and the p- value de- 
creases correspondingly. The Table compares fit 
probabilities for two different assumptions on the 



^Our fits are relatively insensitive to Aq, so we do not 
display figures for this parameter. 



Higgs boson mass measurements ~ 119, 125 GeV, 
see above, and with the option of dropping the 
((7—2)^ constraint in the latter case0. The combi- 
nation of the increase in a-nd in the increase in 
the number of d.o.f. leads to a substantially lower 
p- value after the inclusion of Eq. ([Ij, if (5 — 2)^ 
is taken into account. On the other hand, a hy- 
pothetical mass measurement at 119 GeV would 
yield an improvement in the fit. For compari- 
son, we also show the parameters for the best- 
fit points. Since the uncertainties are large and 
highly non-Gaussian, we omit them from the Ta- 
ble. 

The restrictions that the hypothetical LHC Alh 
constraint imposes on mi/2 are also visible in 
Fig. [3l where we display the effects of an LHC 
Mh constraint in the (mi/2,tan/3) planes of the 
CMSSM and NUHMl. We see here that an LHC 
Alh constraint enlarges visibly the 68% CL range 

''The fit probabilities are indicative of the current exper- 
imental data preferences for one scenario over another 
but, as discussed in [3], but they do not provide a robust 
confidence-level estimation for the actual choice made by 
Nature. 
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Table 1 

Comparison of the best- fit points found in the CMSSM and NUHMl pre-Higgs and for the two potential 
LHC Higgs mass measurements discussed in the text: Mh — 119 and 125 GeV. In the latter case, we also 
quote results if the {g — 2)^ constraint is dropped. At the best-fit NUHMl points, the common values of 
the soft supersymmetry-breaking contributions to the Higgs squared masses are the following - pre-Higgs: 
-1.2 X 10^ GeV^ with Mh ~ 125 GeV and {g - 2)^; -5.5 x 10^ GeV^, with Mh ~ 125 GeV hut without 
[g - 2)^; -8.6 x lO"^ GcV^, with Mh ^ 119 GeV and {g - 2)^; -1.2 x 10^ GeV^ 



of tan/3 in the NUHMl, whereas the change is 
less pronounced in the CMSSM. 

The results for the (M^,tan/3) planes in the 
CMSSM and the NUHMl are shown in Fig. |2 
We observe a strong increase in the best-fit value 
of Ma in both models, especially in the CMSSM, 
where now Ma — 1600 GeV is preferred. We re- 
emphasize, however, that the likelihood function 
varies relatively slowly in both models, as com- 
pared to the pre-LHC fits. 

We now discuss the CMSSM and NUHMl pre- 
dictions for some of the most interesting super- 
symmetric obscrvablcs for the LHC in light of a 
possible LHC measurement at Mh — 125 GeV. 

The upper panels of Fig. [5] display the one- 
dimensional functions for rUg before and af- 
ter applying the new LHC Mh ~ 125 GeV con- 
straint (dashed and solid lines, respectively, in 
both cases including {g — 2)^). We also show as 
dotted lines the Ax^ functions for a fit includ- 



ing Mh — 125 GeV and dropping [g — 2)^. As 
expected on the basis of Fig. [21 the preferred val- 
ues rrig ~ 4 TeV in the CMSSM are much higher 
than in our pre-LHC fit and what would be pre- 
ferred \i Mh — 119 GeV, and presumably beyond 
the reach of the LHC. On the other hand, in 
the NUHMl rrig ^ 2 TeV is marginally preferred. 
However, in both models the Ax^ function varies 
little over the range (2, 4) TeV. Similar features 
are found for niq^^ , as shown in the lower pan- 
els of Fig. m In both models, the regions of rUg 
and TOq^ with Ax^ ^ 1 start at masses around 
1.5 TeV, leaving a large range accessible to the 
SUSY searches at the LHC. In the case of the 
lighter stau mass m^^ for Mh — 125 GeV shown 
in Fig. [6l we again see preferred masses larger 
than in the pre-Higgs fit, with favoured values 
extending up to mf^ ~ 1 TeV. 

We now turn to the predictions of our fits for 
BR(i?s — > shown in Fig. [T] This observ- 
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Figure 3. The (mi/2,tan/3) planes in the CMSSM (left) and the NUHMl (right), for Mh ~ 125 GeV. 
The notations and significations of the contours are the same as in Fig. [H 
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Figure 4. The (A/a, tan /3) planes in the CMSSM (left) and the NUHMl (right), for Mh 
notations and significations of the contours are the same as in Fig. [H 



125 GcV. The 



able is not very sensitive directly to Mh, and the 
indirect sensitivity via mi/2 is not very strong, 
though smaller values of mi/2 do lead to larger 
values of BR(_Bs — > in general. As seen in 

Fig. [3 imposing the putative LHC constraint 
indeed has little effect on BK{Bs fi^ij,~). We 
recall that the best-fit values in the CMSSM and 
NUHMl are both slightly larger than in the SM, 
and enhancements of up to 0(30 — 40%) with re- 
spect to the SM prediction could be detected at 
the LHC at the 3 a level. 



Finally, in Fig. |8] we show results for the pre- 
ferred regions in the (771^^0,17^^) plane. As seen 
in Fig. [HI the fact that larger values of mi/2 
and hence m-^o arc favoured by the larger val- 
ues of Mh implies that at the 68% CL the pre- 
ferred range of a^^ is significantly lower when 
Mh — 125 GeV, when compared to our previous 
best fit with Mh = 119 GeV, rendering direct de- 
tection of dark matter significantly more difficult. 
Again, this effect on m^o is more pronounced in 
the CMSSM, whereas in the NUHMl the value of 
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Figure 5. The one- dimensional Ax^ functions for rrig (upper) and rriq^^ (lower) in the CMSSM (left) and 
the NUHMl (right). The solid lines are for fits assuming Mh — 125 GeV and including [g — 2)^, and the 
dotted lines for fits with ~ 125 GeV but without {g — 2)^. The dashed lines show the results for fits 
without Mh — 125 GcV hut including [g — 2)^ 



m^o for the best-fit point changes only slightly. 

Results dropping the {g — 2)^ constraint 

We have restricted our attention so far to Mh — 
125 GeV assuming the {g — 2)^ constraint. How- 
ever, this value of Mh corresponds approximately 
to our best-fit points in [3] when the {g — 2)^ 
constraint is dropped 0. Accordingly, we now 
consider an the same measurement as given in 



^We recall that it was shown in [3] that the 
CMSSM/NUHMl interpretation of (g - 2)^ is in some 
tension with the LHC constraints on events withl^T- 



Eq. ([T]), but with {g — 2)^ dropped from the fit0. 
In the following plots we show results for fits omit- 
ting {g — 2),,, pre-Higgs (dotted) and post-Higgs 
(solid). 

We see in Fig. |9] that the regions of the 
(mo,mi/2) planes in the CMSSM and NUHMl 
that are favoured at the 68% CL are concen- 
trated at large values if the {g — 2)^ constraint 
is dropped. This reflects the relative harmony 

^There are small differences between the pre-Higgs 68 and 
95% CL contours presented here and the corresponding 
contours in [3] , which provide a measure of the uncertain- 
ties in the interpretation of the MCMC data generated for 
our analysis. 
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Figure 6. The one- dimensional Ax^ functions for the lighter scalar tau mass m^^ in the CMSSM (left) 
and the NUHMl (right), for Mfi ~ 125 GcV. The notations and significations of the lines are the same 
as in Fig. [31 
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Figure 7. The one- dimensional Ax^ functions for IiIl{Bs — > ) in the CMSSM (left) and the NUHMl 

(right), for Mh — 125 GeV. The notations and significations of the lines are the same as in Fig. O 



between the LHC constraints and the hypo- 
thetical Mh — 125 GeV measurement if {g — 2)^j 
is omitted. The inclusion of Eq. ([T|) substantially 
sharpens the prediction at the 68% CL, whereas 
it is less pronounced for the 95% CL contours. 

As we see in Fig. [TUl the concentration at rela- 
tively large mi/2 is reflected in a correlated pref- 
erence for large values of tan (3. Furthermore, as 
seen in Fig. [TTl the corresponding preferred range 
of Ma is also concentrated at relatively large 
masses. Again the inclusion of Eq. ([Ij consid- 



erably sharpens the preferred parameter ranges. 

Looking back now at the one-dimensional Ax^ 
functions for the fits without (g — 2)^ that are 
shown as dotted lines in Figs. [5] and IH we see 
that the preference for large values of (mo, mi/2) 
carries over into relatively large values of m^, mgj, 
and mf J . In particular, the {g — 2)^-less scenar- 
ios offer somewhat gloomy prospects for sparticle 
detection at the LHC. On the other hand, as seen 
in Figs. [3 there is little change in the prediction 
for BR{Bs if (ff — 2)^ is omitted. 
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Figure 8. The (m^o,cr^i) planes in the CMSSM (left) and the NUHMl (right), for Mh ~ 125 GeV. The 
notations and significations of the contours are the same as in Fig. [B 
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Figure 9. The {mo,mi/2) planes in the CMSSM (left) and the NUHMl (right), for Mh ~ 125 GcV, hut 
dropping the {g — 2)^ constraint. Dotted lines show the contours found previously in dropping the 
{g — 2)fj, but without this Mh constraint. Here the open green stars denote the pre-Higgs best-fit points 
(also dropping {g — 2)^i), whereas the solid green stars indicate the new best- fit points. These best- fit points 
are essentially coincident in the NUHMl case. 



Turning finally to the predictions for if [g — 
2)^ is omitted, shown in Fig. [121 wc see that the 
relatively large values of mi/2 seen in Fig. [9] are 
reflected in relatively large values of m^o^, which 
correspond in turn to relatively low values of (7^^. 
The inclusion of Eq. ([T} again strongly reduces 



the preferred parameter ranges. 

An alternative interpretation of a Higgs sig- 
nal around Mh — 125 GeV would be that while 
the MSSM might still be realized, it is not the 
CMSSM nor the NUHMl that describes Nature 
correctly, but another version of the MSSM. In 
this case, the prospects for sparticle detection at 
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Figure 11. The {MA,ta.n/3) planes in the CMSSM (left) and the NUHMl (right), for M,, ~ 125 GcV, 
but dropping the (g — 2)^ constraint. The notations and significations of the contours are the same as in 
FigM 



the LHC and dark matter detection might weh be 
more cheerful than in the (g — 2)^-less CMSSM 
and NUHMl scenarios discussed here. However, 
the exploration of such possible alternative mod- 
els lies beyond the scope of our analysis. 

What if Mft = 119 GeV? 

We have restricted our attention so far to Mh ~ 



125 GeV, which corresponds to the excess seen 
in both CMS and ATLAS. We now consider an 
alternative potential LHC measurement Mh = 
119±1 GeV, which corresponds to the CMS ZZ* 
signal and our earlier predictions including the 
[g — 2)^ constraint, again allowing for a theoret- 
ical error ±1.5 GeV in the calculation of Mh for 
any given set of CMSSM or NUHMl parameters. 
The (mo,TOx/2) planes shown in Fig. [T3] for 
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Figure 12. The {m^o,a^^) planes in the CMSSM (left) and the NUHMl (right), for Mh ~ 125 GeV, hut 
dropping the {g — 2)^ constraint. The notations and significations of the contours are the same as in 
FigM 



the CMSSM (left) and NUHMl (right), the pre- 
ferred regions are shown at the 68% CL (red) 
and 95% CL, with the sohd (dotted) hnes include 
(omit) the assumed LHC Higgs constraint. Since 
this assumed LHC value of Mh coincides with the 
previous best-fit values in both the CMSSM and 
NUHMl, the best-fit points in these models (in- 
dicated by the green stars in Fig. [T5)) are nearly 
unaffected by the imposition of the putative LHC 
constraint. The effect of the hypothetical mea- 
surement restricting the range in TOi/2 is indeed 
seen in both panels of Fig. [131 though for the 68% 
CL contour (shown in red) it is much more pro- 
nounced for the CMSSM than for the NUHMl, 
whereas for the 95% CL contour (shown in blue) it 
is more significant for the NUHMl. This reflects 
the fact in the NUHMl the global Ax^ function 
found in [3] rose quite steeply in the neighbour- 
hood of the best-fit point, resulting in a relatively 
tight 68% CL contour, whereas the rise of away 
from the best-fit point in the CMSSM was more 
gradual. This led previously to a larger 68% CL 
contour and a broader range of Mh at the 68% 
CL, which is now more affected by an assumed 
LHC Mh constraint. On the other hand, the 95% 
CL contour in the NUHMl extended previously 
to larger values of mi/2 than in the CMSSM, and 
these values arc particularly susceptible to the 



LHC Mh constraint. 

Although we add another constraint (as dis- 
cussed above), the total at the best-fit points 
do not change 0. For this reason, the p- values for 
the CMSSM and NUHMl would increase for a 
hypothetical measurement Mh ~ 119 GeV, cor- 
responding formally to better overall fits to the 
larger data set, as seen in Table [TJ 

As one might expect, such an LHC A4h con- 
straint would reduce considerably the 68% CL 
range of tan/3 in the CMSSM. This is because, 
for mi/2 close to the best-fit value, ^ 700 to 
800 GeV, fixing the Higgs mass at 119 GeV dis- 
favours low values of tan /?, which yield low val- 
ues of Mh- This effect is not important in the 
NUHMl, where the range of tan/3 was already 
smaller before imposing the Mh constraint. We 
also note that Ma is restricted to somewhat 
smaller values when the hypothetical LHC con- 
straint on Mh is included. Furthermore, as ex- 
pected, the values of nrig at the minima of the 
functions are not affected and there is lit- 
tle change in for rUg between 2 and 3 TeV. 
(The corresponding plots are not shown.) How- 
ever, there are significant effects at both lower 

^They would change only slightly if the Higgs mass were 
assumed to differ by < 1 GeV from that obtained at the 
best-fit point. 
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Figure 13. The (^0,^1/2) planes in the CMSSM (left) and the NUHMl (right) assuming a hypothetical 
measurement of Mh — 119 GeV. The notations and significations of the contours are the same as in 
Fig.\E 
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Figure 14. The {m^o,a^^) planes in the CMSSM (left) and the NUHMl (right), for Mu 
notations and significations of the contours are the same as in Fig. \13\ . 



119 GcV. The 



and higher values of m^. 



values of 3 TeV are disfavoured. 



In particular, large 
The 

prospects for discovering gluinos at the LHC in 
the near future would remain uncertain in both 
the CMSSM and NUHMl. An LHC measure- 
ment of Mh — 119 GeV would disfavour large 
squark masses, but the 95% CL range would still 
extend to m^^ ~ 4 TeV in the CMSSM and 
~ 2 TeV in the NUHMl. The preferred value 



of rrif, - 300 GeV in both the CMSSM and 
NUHMl both with and without the hypothetical 
LHC Mh measurement, with large masses again 
becoming somewhat more disfavoured. 

Finally, in both the CMSSM and the NUHMl 
there is little impact on the 95% CL regions nor 
on the 68% CL region in the NUHMl in the 
(wm-o I CTp') plane. The only substantial change, 
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as can be seen in Fig. [TH appears in the 68% 
CL region of the CMSSM, where now values of 
TO^o > 700 GeV and cr^i < lO^^'^cm-^ are dis- 
favoured after the inclusion of a Higgs-boson mass 
measurement at 119 GeV. 



4. Summary and Conclusions 

The ATLAS and CMS searches for the Higgs 
boson have already excluded a very large range 
of masses, with the only remaining windows for 
a SM-like Higgs boson being in the ranges Mh G 
(115.5, 127) GeV or > 600 GeV [23l[24]. The lat- 
ter range is disfavoured by precision electroweak 
data, so attention naturally focuses on the low- 
mass range. It may or not be a coincidence that 
this range includes the range Mh $130 GeV ac- 
cessible in simple supersymmetric models such as 
the CMSSM and NUHMl. Within this range, 
our previous global fits of these models includ- 
ing \g - 2)^ predicted Mh 119 GeV if the 
{g — 2)^ constraint was included in the fit, and 
Mh ^ 126 GeV if {g - 2)^ was omitted [3]. The 
latest ATLAS and CMS results display an inter- 
esting fiuctuation at Mh ^ 125 GeV (i.e. close to 
the latter result from [3]) and we have combined 
a hypothetical measurement of Mh ~ 125 GeV 
with the global likelihood functions obtained in 
our previous fits p]- 

As we have shown in this paper, this combi- 
nation refines our previous predictions for the 
CMSSM and NUHMl model parameters within 
global fits incorporating [g — 2)^. In particular, 
the combination prefers a range of larger values of 
mi/2, resulting in larger values of rrig and other 
sparticle masses being preferred, restricting the 
prospects for discovering supersymmetry at the 
LHC within these models. The predictions for 
CTp^ are pushed to higher masses and lower cross 
sections, particularly in the CMSSM. There are 
also smaller changes in the predictions for other 
observables such as BR(i3s — > . 

We have also shown the analogous CMSSM and 
NUHMl fit results for a hypothetical measure- 
ment of Mh — 125 GeV if the (g — 2)^ constraint 
is omitted. In this case we find a stronger pref- 
erence for larger values of {mQ,mi/2), and cor- 
respondingly larger values of tan/3 and Ma, as 



well as larger values of mg,mg^, potentially ly- 
ing beyond the reach of the LHC. We have also 
commented on the potential implications of a hy- 
pothetical Higgs discovery at Mh ~ 119 GcV. 

Time will soon tell where the LHC experiments 
are indeed discovering the Higgs boson. How- 
ever, we have shown that Mh = 125 GeV is a 
possibihty within the CMSSM and NUHMl, al- 
though it lies at the upper range of what is pos- 
sible within the CMSSM or NUHMl, and might 
suggest reduced prospects for discovering these 
particular models of supersymmetry at the LHC. 
Alternatively, it could well be that one should 
look beyond the frameworks of the models dis- 
cussed here. 

Note Added 

After acceptance of this paper for publication, 
we became aware of issues in the implementation 
of the FeynHiggs code and in the cold dark mat- 
ter density calculation, which required extra sam- 
pling and reprocessing of the NUHMl parameter 
space. We arc grateful to Nazila Mahmoudi and 
Azar Mustafayev for discussions on dark matter 
density calculations. 

Acknowledgements 

The work of O.B., K.J.D., J.E., J.M. and 
K.A.O. is supported in part by the London Cen- 
tre for Terauniverse Studies (LCTS), using fund- 
ing from the European Research Council via the 
Advanced Investigator Grant 267352. The work 
of S.H. is supported in part by CICYT (grant 
FPA 2010~22163-C02-01) and by the Spanish 
MICINN's Consolider-Ingenio 2010 Program un- 
der grant MultiDark CSD2009-00064. The work 
of K.A.O. is supported in part by DOE grant 
DE-FG02-94ER-40823 at the University of Min- 
nesota. 

REFERENCES 

1. J. R. Ellis, S. Heinemeyer, K. A. Olive and 
G. Weiglcin, Phys. Let t. B 515 (2001) 348 
[arXiv:hep-ph/0105061| . 

2. S. Ambrosanio, A. Dedes, S. Heinemeyer, 
S. Su and G. Weiglein, Nucl. Phys. B 624 



(2002) 3 [arXiv:hep-ph/0106255| . 



16 



6. 



7 



O. Buchmucllcr, et al., Eur. Phys. J. C 72, 
1878 (2012) I ar Xiv: 11 10.35681 [hep-ph]]. 
J. R. Ellis, G. Ridolfi and F. Zwirner, Phys. 
Lett. B 257 (1991) 83; Phys. Lett. B 262 
(1991) 477; Yasuhiro Okada, Masahiro Yam- 
aguchi and Tsutomu Yanagida, Phys. Lett. 
B262, 54, 1991; Prog. Theor. Phys. 85 (1991) 
1; A. Yamada, Phys. Lett. B 263 (1991) 
233; Howard E. Haber and Ralf Hempfling, 
Phys. Rev. Lett. 66 (1991) 1815; M. Drees and 
M. M. Nojiri, Phys. Rev. D 45 (1992) 2482; 
P. H. Chankowski, S. Pokorski and J. Rosiek, 
Phys. Lett. B 274 (1992) 191; Phys. Lett. B 
286 (1992) 307. 

G. Degrassi, S. Heinemeyer, W. Hollik, 
P. Slavich and G. Wciglcin, Eur. Phys. J. C 28 
(2003) 133 |arXiv:hcp-ph/0212020 ; S. Heine- 
meyer, W. Hollik and G. Weiglein, Eur. Phys. 
J. C 9 (1999) 343 arXiv:hep-ph/9812472j; 
S. Heinemeyer, W. Hollik and G. Wei- 
glein, Comput. Phys. Commun. 124 
(2000) 76 arXiv:hep-ph/9812320'; 
M. Fra nk et ai, JHEP 702 (2007) 
047 [arXiv:hep-ph/0611326| ; See 
|htt p : / / www ■ f eynhiggs . del . 
M. Carcna, P. Draper, S. Heinemeyer, T. Liu, 



C. E. M. Wagner and G. Weiglein, Phys. Rev. 
D 83 (2011) 055007 (arXiv: 101 1. 5304 [hep- 
ph]]. 

O. Buchmueller et al, Phys. Lett. B 657 
(2007) 87 larXiv:0707.3447' [hep-ph]]. 

8. O. Buchmueller et ai, JHEP 0809 (2008) 117 
|arXiv:0808.4128l [hep-ph]] . 

9. O. Buchmueller et al, Eur. Phys. J. C 64 
(2009) 391 arXiv:09QL5568| [hep-ph]]. 

10. O. Buchmueller et al, Phys. Rev. D 81 (2010) 
035009 arXiv:0912.1036 [hep-ph]]. 
O. Buchmueller et al, Eur. Phys. J. C 71 
(2011) 1583 jarXiv: 101 1.61 181 [hep-ph]]. 
O. Buchmueller et al, Eur. Phys. J. C 71 
(2011) 1634 jarXiv:1102.4585l [hep-ph]]. 
O. Buchmueller et al, Eur. Phys. J. C 71 
(2011) 1722 [arXiv:1106.2529l [hep-ph]]. 

14. For more information and updates, please see 



11 



12 



13 



http : //cern . ch/mastercode/ 



15. For a sampling of other pre-LHC analyses, 
see: E. A. Baltz and P. Gondolo, JHEP 0410 
(2004) 052 |arXiv:hep-ph/0407039|; B. C. Al- 



lanach and C. G. Lester, Phys. Rev. D 73 
(2006) 015013 [arXiv:liep-ph/0507283] ; 
R. R. de Austri, R. Trotta and 
L. Roszkowski, JHEP 0605 (2006) 002 
farXiv:hep-ph/060202g; R. Lafaye, T. Plehn, 
M. Ranch and D. Zerwas, Eur. Phys. J. 
C 54 (2008) 617 ;arXiv:0709.3985' [hep- 
ph]]; S. Heinemeyer, X. Miao, S. Su and 
G. Weiglein, JHEP 0808 (2008) 087 
[arXiv:0805.2359' [hep-ph]]; R. Trotta, 
F. Feroz, M. P. Hobson, L. Roszkowski and 
R. Ruiz de Austri, JHEP 0812 (2008) 024 
[arXiv:0809.3792 [hep-ph]]; P. Bechtle, K. De- 
sch, M. Uhlenbrock and P. Wienemann, Eur. 
Phys. J. C 66 (2010) 215 [arXiv:0907.2589l 
[hep-ph]]. 

16. For a sampling of other post-LHC analy- 
ses, see: D. Feldman, K. Frecse, P. Nath, 
B. D. Nelson and G. Peim, Phys. Rev. 
D 84, 015007 (2011) |arXiv:1102. 25481 [hep- 
ph]]; B. C. All anach, Phys. R ev. D 83, 
095019 (2011) gar Xiv:1102.3T4 9' [hep-ph]]; 
S. Scopel, S. Choi, N. Fornengo and 
A. Bo ttino, Phys. Re v. D 83, 095016 
(2011) [arXiv:1102.4033 [hep-ph]]; P. Bech- 
tle et al, Phys. Rev. D 84, 011701 (2011) 
|arXiv:1102.4693l [hep-ph]]. B. C. Allanach, 
T. J. Khoo, C. G. Lester an d S. L. WiUiams, 
JHEP 1106, 035 (2011) l arXiv: 1103. 09691 
[hep-ph]]; S. Akula, N. Chen, D. Feldman, 
M. Liu, Z. Liu, P. Nath and G. Peim, Phys. 
Lett. B 699, 377 (2011) larXiv: 1103. 11971 
[hep-ph]]; M. J. Dolan, D. Grellscheid, 
J. Jaeckel, V. V. Khoze and P. Richard- 
son, JHEP 1106, 095 (2011) [arXiv:llQ4,_6585l 
[hep-ph]]; S. Akula, D. Feldman, Z. Liu, 
P. Nath and G. Peim, Mod. Phys. Lett. A 
26, 1521 (2011) |arXiv:1103.5 '06n [hep-ph]]; 
M. Farina, M. Kadastik, D. Pappadopulo, 
J. Pata, M. Raidal and A. Strumia, Nucl. 
Phys. B 853, 607 (2011) [arXiv:1104.3572l 
[hep-ph]]; S. Pr ofumo, Phys. R ev. D 84, 
015008 (2011) |arXiv:1105.5T62l [hep-ph]]; 
T. Li, J. A. Maxin, D. V. Nanopoulos 
and J. W. Walker, larXiv:1106.lT65l [hep- 
ph]; N. Bhattacharyya, A. Choudhury and 
A. Datta, Phys. Rev. D 84, 095006 (2011) 
larXiv: 1107. 19971 [hep-ph]] . 



17 



17. A. Fowlie, A. Kalinowski, M. Kazana, 
L. Roszkowski and Y. L. S. Tsai, Phys. Rev. D 
85, 075012 (2012) |arXiv:1111.6 098 [hep-ph]]. 

18. R. Baratc et al. [ALEPH, DELPHI, L3, 
OPAL Collaborations and LEP Working 
Group for Higgs boson searches], Phys. Lett. 
B 565 (2003) 61 [arXiv:hep-ex/0306033'. 

19. S. Schacl et al. [ALEPH, DELPHI, L3, OPAL 
Collaborations and LEP Working Group for 
Higgs boson searches], Eur. Phys. J. C 47 
(2006) 547 [arXiv:hep-ex/0602042 . 

20. |http : //tevnphwg . f nal . gov/ and references 
therein. 

21. ATLAS Collaboration, 

https : //atlas . web . cern . ch/Atlas/GROUPS/ 

PHYSICS/CONFNOTES/ATLAS-CONF-2011-132/ 

ATLAS-C0NF-2011-132.pdf. 

22. CMS Collaboration, 



[http : / / cdsweb . cern . ch/ record/ 1378096/ 
f iles/HIG-ll-020-pas.pdf . 

23. CMS Collaboration, 

https : // cdsweb . cern. ch/ ecord/ 1399607/ 
f iles/HIG-ll-023-pas.pdf ; 
ATLAS Collaboration, 

https : / / cdsweb . cern. ch/record/1399599 

?ln=enATLAS-C0NF-2011-157; 

L. Rolandi, on behalf of ATLAS, CMS 

and the LHC Higgs Combination Group, 

|http : //indico . in2p3 . f r/ getFile . py/ 

access?contribId=72&sessionId=19&resId 

=0&materialId=slides&conf Id=6004. 

24. F. Gianotti for the ATLAS Collabora- 
tion, G. Tonelli for the CMS Collaboration, 
[http : / / indico . cern . ch/conf erenceDisplay 

.py?conf Id=164890. 

25. See, for example, J. Ellis, J. R. Espinosa, 
G. F. Giudice, A. Hocckcr, A. Riotto, Phys. 
Lett. B 679 (2009) 369 and references therein. 

26. J. Ellis, D. A. Ross, Phys. Lett. B 506 (2001) 
331. 

27. M. Drees and M. M. N ojiri , Phys. Rev . 
D 47 (1993) 376 |arXiv:hep-ph/9207234] ; 
G. L. Kane, C. F. Kolda, L. Roszkowski 
and J. D. Wells, Phys. Rev. D 49 (1994) 
6173 |arXiv:hep-ph/9312272i ; H. Baer 
and M. Brhlik, Phys. Rev. D 53 (1996) 
597 |arXiv:hcp-ph/950832l; Phys. Rev. 
D 57 (1998) 567 |arXiv:hep-ph/9706509l ; 



J. R. Elhs, T. Falk, K. A. Olive and 
M. Schmitt, Phys. Lett B 388 (1996) 
97 |arXiv:hep-ph/9607292|; Phys. Lett. B 
413 (1997) 355 "arXiv:hcp-ph/970544^; 
J. R. Ellis, T. Falk, G. Ganis, K. A. Olive 
and M. Schmitt, Phys. Rev. D 58 (1998) 
095002 'arXiv:hep-ph/9801445l; V. D. Barger 
and C. Kao, Phys. Rev. D 57 (1998) 3131 



arXiv:hep-ph/9704403| ; J. R. Elhs, T. Falk, 
G. Ganis and K. A. Olive, Phys. Rev. D 
62 (2000) 075010 arXiv:hep-ph/0004169 ; 
J. R. Ellis, T. Falk, G. Ganis, K. A. Olive 
and M. Srcdnicki, Phys. Lett. B 510 (2001) 
236 arXivdiepjj3h/0102 098| ; V. D. Barger 
and C. Kao, PhysrLet t. B 518 (2001) 117 
[arXiv:hep-ph/0106189'; L. Roszkowski, 
R. Ruiz dc Austri and T. Nihei, JHEP 
0108 (2001) 024 [arXiv:hep-ph/0106334] ; 
A. Djouadi, M. Drees and J. L. Kneur, JHEP 
0108 (2001) 055 arXiv:hep-ph/0 107316 ; 
U. Chattopadhyay, A. Corsetti and 
P. Nath, Phys. Rev. D 66 (2002) 035003 
arXiv:hcp- ph/0201001 ; J. R. EUis, 
K. A. Olive and Y. Santoso, New Jour. Phys. 
4 (2002) 32 |arXiv:hep-ph/0202110| ; H. Baer, 
C. Balazs, A. Belyaev, J. K. Mizukoshi, 
X. Ta ta and Y. W ang, JHEP 0207 (2002) 
050 |arXiv:hep-ph/0205325; R. Arnowitt 
and B. Dutta, arXiv:hcp-ph/0211417^ 
J. R. EUis, K. A. Olive, Y. Santoso and 
V. C. Spanos, Phys. Lett. B 565 (2003) 
176 arXiv:hep-ph/0303043| ; H. Baer 
and C. Balazs, JCAP 0305, 006 (2003) 
arXiv:hep-ph/0303114"; A. B. Lahanas 
and D. V. Nanopoulos, Phys. Lett. B 
568, 55 (2003) |arXi^hep-ph/0303130] ; 
U. Chattopadhyay, A. Corsetti and 
P. Nat h, Phys. Rev. D 68, 035005 
(2003) [arXiv:hep-ph/0303201| ; C. Munoz, 
Int. J . Mod. Phys. A 19, 3093 (2004) 

Arnowitt, 



R. 



|arXiv:hep-ph/03 934 6 ; 
B. Dutta and B. Hu, |arXiv:hep-ph/03l0103 
J. Ehis and K. A. Olive, in "Particle 
Dark Matter", ed. G. Bertone, ppl42-163 
[arXiv:1001.3651l [astro-ph.CO]] . 
28. J. R. Elhs, D. V. Nanopoulos, K. A. Olive and 
Y. Santoso, Phys. Lett. B 633, 583 (2006) 
[hep-ph/050933ll. 



18 



29. H. Bacr, A. Mustafaycv, S. Profumo, 
A. Belyaev and X. Tata, Phys. Rev. D 
71 (2005) 095008 |j arXiv:hep-ph/0412059| ; 
H. Baer, A. Mustafayev, S. Profumo, 
A. Belyaev and X. T ata, JHEP 0507 (2005) 
065, |hep-ph/0504001 ; J. R. Ellis, K. A. Olive 

(2008) 



78 



and P. Sandick, Phys. Rev. D 
075012 ^arXiv:0805.2343 [hcp-ph]]. 

30. E. Aprile et al. [XENONIOO Collabora- 
tion], Phys. Rev. Lett. 107, 131302 (2011) 
[arXiv:1104.2549 [astro-ph.CO]]. 

31. S. Heinemeyer et al, J HEP 0608 (2006) 
052 farXiv:hep-ph/0604147 ; S. Heinemeyer, 
W. Hollik, A. M. Weber a nd G. Weiglein, 
JHEP 0804 (2008) 039 |arXiv: 0710 .29721 
[hep-ph]]. 

32. B. C. AUanach, Comput. Phys. Commun. 143 
(2002) 305 [arXiv:hep-ph/0104145j . 

33. G. Isidori and P. Paradisi, Phys. Lett. 
B 639 (2006) 499 arXiv:hep-ph/06 05012 ; 
G. Isidori, F. Mescia, P. Paradisi and 
D. Temes, Phys. Rev. D 75 (2007) 
115019 arXiv:hep-ph/0703035] , and refer- 
ences therein. 

34. F. Mahmoudi, Comput. Phys. Commun. 
178 (2008) 745 |arXiv:0710.2 067 [hep-ph]]; 
Comput. Phys. C ommun. 180 (2009) 1579 
|arXiv:0808.3144' [hep-ph]]; D. Eriksson, 
F. Mahmoudi and O. Stal, JHEP 0811 (2008) 
035 [a rXiv:0808.3551l [hep-ph]] : A. Arbey and 
F. Mahmoudi, Comput. Phys. Commun. 181 
(2010) 1277 jarXiv:0906.0 369 [hep-ph]]. 

35. G. Belanger, F. Boudjema, A. Pukhov and 
A. Semenov, Comput. Phys. Commun. 
176 (2007) 367 |arXiv:hep-ph/06 07059 ; 
Comput. Phys. Commun. 149 (2002) 
103 |arXiv:hep-ph/0112278| ; Com- 
put. Phys. Commun. 174 (2006) 577 



arXiv:hep-ph/0405253| . 



36. Information about this code is available from 
K. A. Olive: it contains important contri- 
butions from T. Falk, A. Ferstl, G. Ganis, 
A. Mustafayev, J. McDonald, K. A. Olive, 
P. Sandick, Y. Santoso and M. Srednicki. 

37. P. Skands et al, J HEP 0407 (2004) 036 
farXiv:hep-ph/0311123 ; B. Allanach et al, 
Comput. Phys. Commun. 180 (2009) 8 



38. G. Aad et al. [ATLAS Collaboration], Phys. 
Lett. B 710, 67 (2012) arXiv: 1109.6572 [hep- 
ex]]. 

39. S. Chatrchyan et al. [CMS Cohabora- 
tion], Phys. Rev. Lett. 107, 221804 (2011) 
jarXiv:1109. 23521 [hep-ex]]. 

40. S. Chatrchyan et al. [CMS CoUabora- 
tion], Phys. Rev. Lett. 107, 191802 (2011) 
| arXiv:1107.5834l [hep-ex]]. 

41. R. Aaij et al. [LHCb Collaboration], Phys. 
Lett. B 699 (2011) 330 arXiv: 1103.246^ 
[hep-ex]]; Phys. Lett. B 708, 55 (2012) 
[arXiv:1112.1600' [hep-ex]]. 

42. T. Aahonen et al. [CDF Collaboration], Phys. 
Rev. Lett. 107, 239903 (2011) [Phys. Rev. 
Lett. 107, 191801 (2011)] |arXiv:1107.2304l 
[hep-ex]]. 

43. S. S. Abdus Salam, et al., Eur . Phys. J. C 71, 
1835 (2011) [arXiv:1109.3859' [hep-ph]]. 

44. P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, 
M. Schelkc and E . A. Baltz, JCA P 0407 
(2004) 008 arXiv:astro-ph/0406204 ; 

[http : //www .physto . se/ ^ edsjo/darksusy/| 

45. See the last reference in [34] . 



jarXiv:0801. 00451 [hep-ph]] 



